Fixed-iV Superconductivity: 
The Exact Crossover from the Bulk to the Few-Electron Limit 
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We use two truly canonical approaches to describe superconductivity in ultrasmall metallic 
grains: (a) a variational fixed- N projected BCS-like theory and (b) the exact solution of the model 
Hamiltonian developed by Richardson in context with Nuclear Physics. Thereby we obtain a de- 
scription of the entire crossover from the bulk BCS regime (mean level spacing d <C bulk gap A) 
to the "fluctuation-dominated" few-electron regime (d > A). A wave-function analysis shows in 
detail how the BCS limit is recovered and how for d > A pairing correlations become delocalized 
in energy space. 



I. INTRODUCTION 



When a system of (correlated) electrons is sufficiently 
small, the electronic spectrum becomes discrete. Such 
a discrete spectrum was directly measured for the first 
time by Ralph, Black and Tinkham (RBT) for ul- 

trasmall Al grains. This allowed them to study the na- 
ture of pairing correlations in a small superconductor in 
unprecedented detail. These experiments gave new actu- 
ality to an old and fundamental question: What is the 
lower size limit for superconductivity? Besides Ander- 
son's prediction || that superconductivity breaks down 
once the single-particle mean level spacing d becomes 
larger than the (bulk) superconducting gap A, the fi- 
nite size of a superconducting grain also manifests itself 
in its large charging energy, which effectively suppresses 
electron number fluctuation on the grain. Hence super- 
conductivity on small grains has to be formulated in a 
manifestly canonical way. 

After briefly introducing the experiments and a toy 
model which captures their essential features, we show 
how the entire crossover between the bulk BCS-like 
regime and the few-electron regime can be described the- 
oretically by fixing the electron number on the grain and 
using either a projected BCS approach or, even better, 
a long-forgotten exact solution of the reduced BCS model 
Hamiltonian. Both the projected BCS approach and the 
exact solution enable us (i) to significantly improve pre- 
vious g.c. upper bounds on ground state energies, [^-0, 
in the latter case by giving the exact result; (ii) to find 
in the crossover regime a remnant of the "break-down 
of superconductivity" obtained in g.c. studies, at which 
the condensation energy changes from being extensive to 
practically intensive; and (iii) to study this change by 
an explicit wave-function analysis, which shows in detail 
how the BCS limit is recovered for d« A, and how for 
d A pairing correlations delocalize in energy space. 



Spectroscopic Gap in an Ultrasmall 
Superconducting Grain 



In RBT's experiments, an ultrasmall grain was used 
as central island in a single-electron transistor: it was 
connected via tunnel barriers to external leads and ca- 
pacitively coupled to a gate, and its electronic spectrum 
determined by measuring the tunnel current through the 
grain as a function of transport voltage at a fixed tem- 
perature of 50mK. For a typical grain the radius was 
r w 5nm, and the crude order-of-magnitude free-electron 
estimate d = 2tt 2 Ti 2 / '(rnkpYol) for the mean level spacing 
near sf yields d ~ 0.5meV. The grain's charging energy 
was about Ec = e 2 /2Ctotai = 50meV and therefore much 
larger than all other energy scales, such as the Aluminum 
bulk gap (A ~ 0.4meV), typical values of the transport 
voltage (V ~ 2mV) and the temperature. 

The most remarkable feature of the experimental tun- 
neling spectrum, shown in Fig. |l|, is the presence of a 
clear spectroscopic gap for the grain with even electron 
number and its absence for the odd grain. This reveals 
the presence of pairing correlations: in even grains, all ex- 
cited states involve at least two BCS quasi-particles and 
hence lie significantly above the ground state, whereas 
odd grains always have at least one quasi-particle and 
excitations need not overcome an extra gap. 

The charging energy, being the largest energy scale of 
the system, strongly suppresses particle number fluctua- 
tions on the grain and hence the discrete energies mea- 
sured in RBT's experiments essentially correspond to the 
eigenspectrum of a grain with fixed electron number N. 
We therefore consider below an ultrasmall grain com- 
pletely isolated from the rest of the world, e.g. by in- 
finitely thick oxide barriers. Our main goal will be to 
elucidate, within a canonical framework, the nature of 
the pairing correlations in an ultrasmall grain at T = 0. 
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II. MODELLING AND FIRST APPROXIMATION 

To investigate the influence of pairing correlations, on 
the excitation spectrum of an ultrasmall grain, we model 
the grain by a reduced BCS-Hamiltonian. It has been 
used before to describe small superconducting grains 
and was phenomenologically successful for d < A 
but probably is unrealistically simple for d 3> A, 
for which it should rather be viewed as toy model: 
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N-1 



JV-1 
j,3'=0 



(1) 



The c'j± create electrons in free time-reversed single- 
particle- in-a-box states \j, ±), with discrete, uniformly 
spaced, doubly degenerate eigenenergies Sj = jd + Sq. 
The interaction scatters only time-reversed pairs of elec- 
trons within ujd of Sf- Its dimensionless strength A is 
related to the two material parameters A and ujd via 
the bulk gap equation sinh(l/A) = cj^/A. We chose 
A = 0.224, close to that of Al §. The level spacing d 
determines the number N = 2u}n/d of levels, taken sym- 
metrically around e f , within the cutoff; electrons outside 
the cutoff remain unaffected by the interaction and are 
thus neglected throughout. 



A. Grand-Canonical BCS Approach 

The most direct and easiest theoretical approach ||] 
to describing an ultrasmall grain simply uses the well- 
known grand-canonical (g.c.) variational BCS ansatz for 
the ground state of an even or odd grain (subscript p = 
or 1, respectively): 
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(2) 



(«? + «? = !) 



Vj and Uj are the amplitudes that level j is doubly oc- 
cupied or empty, respectively. Note that on an odd 
grain one electron necessarily is unpaired; to minimize 
its kinetic energy, it is put at the Fermi energy (£j odd — 
Ef)- Minimizing the energy expectation value £^ c = 

(BCS|//|BCS) (p — 0, 1) with respect to Uj and Vj yields 
the even and odd "gap equations" (at T = 0) : 



•r'E 



•A?(d) 



(3) 



These are solved for the even and odd pairing parameters 
Ao and Ai as a function of level-spacing d. Note that 
Jodd is excluded from the odd sum. 



As predicted in 1959 by Anderson, || it turned out 
II that above a critical level spacing the gap equation 
ceases to have a non-trivial solution: when the sam- 
ple becomes too small, superconductivity breaks down. 
More surprising was the finding |||| that the break- 
down is parity- dependent the odd Ai vanishes already 
at a much smaller level spacing (say df c ) than the even 
Ao (say d^ c ). This is reflected in the condensation 
energy E p = £ p — (F p \H\F p ) , which is measured rela- 
tive to the energy of the respective uncorrelated Fermi 
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Fig. |2|(a) below shows that the critical level spacing above 
which the g.c. results for E p c reduce to zero is parity- 
dependent. 



B. Fluctuations 



While the g.c. BCS ansatz confirms Anderson's 
prediciton, it is, however, clear that the approach has 
two problems: (a) it inherently contains fluctuations of 
the particle number N which are not present in the actual 
grain system due to its large charging energy; and (b) it 
fails to describe superconducting fluctuations for large 
level spacings d ^> A, where it trivially yields A = 0. 

For bulk systems, the fluctuations in TV are negligibly 
small (5N ~ yN -C N), and a more rigorous fixed- TV 
treatment would only correct the BCS ground state en- 
ergy per electron H|l(| by order 1/iV, which vanishes in 
the thermodynamic limit d — > 0. For ultrasmall systems, 
however, precisely such corrections are important and we 
will have to incorporate them. 

On the other hand, pairing correlations (the redistri- 
bution of pairs from below to above sp) can lower the 
condensation energy below zero even when A p = 0. In 
this regime, pairing correlations are traditionally called 
"superconducting fluctuations" , which are evidently not 
captured adequately by the g.c. ansatz (||). Matveev and 
Larkin (ML) Q calculated the energy lowering to loga- 
rithmic order for d ^> A (i.e. A p — 0). However, they are 
known |i) to become important already in the crossover 
regime d ~ A. 

To adequately describe an ultrasmall superconducting 
grain we therefore must go beyond standard BCS theory 
by (a) fixing the particle number N and (b) incorporating 
superconducting fluctuations. These ingredients will al- 
low us to describe the full crossover between the bulk sys- 
tem, which is dominated by BCS superconductivity, and 
the few-electron system, which mainly shows supercon- 
ducting fluctuations. In particular we shall study how, at 
the crossing of the two energy scales d and A, the break- 
down of superconductivity predicted by BCS theory is 
softened and how the fluctuations become dominant. 
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III. THE FULL CROSSOVER FROM SMALL N 
TO THE BULK LIMIT 

A. Fixed-TV Projection 

Our first crossover study adapts a method devel- 
oped by Dietrich, Mang and Pradal for shell modells 
of nuclei with pairing interaction to the case of ultrasmall 
metallic grains. This "projected BCS" (PBCS) method 
also is a variational approach, but projects (before varia- 
tion) the trial wave- function onto a fixed electron number 
N = 2riQ + p (p is the parity) : 

|PBCS)o = 

,2tt AT-1 

C d4>e-^ J] KWVj+4-)|Vac>, (4) 

J ° 3=0 

(In the odd case Y\- again is replaced by cj odd nj^j odd -) 
The integral over (f> performs the projection onto the 
fixed electron pair number no, and C is a normaliza- 
tion constant ensuring o(PBCS|PBCS)o = 1. Again, 
the amplitudes Vj and Uj are found by minimizing the 
energy expectation value of the projected wave-function 
f pbcs = p (PBCS|#|PBCS) p . While in the g.c. case the 
wave-function essentially can be described by a single pa- 
rameter A p , in this case the minimization leads to a set 
of 2hq coupled non-linear equations which include projec- 
tion integrals (for details see |ll|). Following Ref. Jl2| , 
we evaluate all integrals numerically (using fast Fourier 
transform routines). 

In the limit d — > at fixed nod, the PBCS theory re- 
duces to the g.c. BCS theory of Ref. || (proving that the 
latter's TV-fluctuations become negligible in this limit): 
The projection integrals can then be approximated by 
their saddle point values [jl2j; at the saddle, the vari- 
ational equations decouple and reduce to the BCS gap 
equation while the saddle point condition fixes the mean 
number of electrons to be 2no- To check the opposite 
limit of d ^> A where no becomes small, the so-called 
fluctuation-dominated regime, we compared the PBCS 
results for E$ with exact results, finding agreement to 
within, say, 6% for no < 10. 1 This shows that super- 
conducting fluctuations are automatically treated ade- 
quately in the PBCS approach. 

The advantages of the PBCS method relative to the 
g.c. one are (a) the similarity of the trial wave-function 



to the BCS wave function, allowing it to capture the bulk 
limit and, on the other hand, (b) the increase of varia- 
tional degrees of freedom from one (A in the g.c. BCS 
theory) to N — 1 (one for each single particle level minus 
1 for fixing TV) due to the projection onto a fixed particle 
number. These additional degrees of freedom allow the 
method to also capture the superconducting fluctuations 
for d ^S> A. Because it works so well for d <C A and 
d » A, one might hope that it acceptably describes the 
crossover regime d ~ A, too. There it really amounts 
to an uncontrolled approximation, whose quality will be 
checked against exact results in the next section. 



B. Richardson's Exact Solution 

Since the publication of |ll[ we became aware of the 
fact that Richardson had shown already in the mid-60's 
jrjj, also in the context of nuclear physics, that the 
ground state of the Hamiltonian (|l|) actually can (for 
non-degenerate Sj) be found exactly by solving a set of 
no coupled algebraic non-linear equations. 2 

Richardson introduces electron pairs bj = Cj—Cj+, 
which have the commutation relation = — 

2bjbj), and exploits the fact that in the Hilbert space of 
non-singly-occupied states, the modified Hamiltonian 

H = J2^ 3 b}b 3 -\dJ2b]b f . (5) 

3 jj' 

is equivalent to the H of ([!]). Solving H exactly would be 
trivial if the b's would represent true bosons. However, 
they actually are "hard-core bosons" instead. Richard- 
son thus expresses the general ground state \G) N of the 
system as 

\G) N = Yl ^•••Jn )fo] 1 ---fo]„jVac), (6) 

where the sum is explictly restricted to exclude double 
occupancy of pair states. The wave-function (p is found 
by solving the many-body Schrodinger equation for ip. 
Richardson showed that the following ansatz works: 

Mi ■ ■ -ho) oc 5> (fi 2f V | ■ (7) 



1 Because this error refers to the correlation energy, it is 
larger than the error (< 1%) cited in jll] for the total ground 
state energy E p . 

2 The numerical implementation of the exact solution is much 
easier than for the PBCS method, since it does not include 
any projections integrals. 
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Here J^v ^ represents the sum over all permutations of 
1, . . . , no, and the parameters Ek are the solution of the 
coupled algebraic equations 



, no 



2no , 



Xd . — ^ i?/ — ^— ' 2e,- — J?j. 
1=1 1 K 7 - = i j K 



k = l...n . (8) 



The total ground state energy is given by £g xact = 

The discussion of Richardson's exact solution in con- 
text with ultrasmall grains will be the subject of a forth- 
coming paper. 



C. Ground State Energies 

Figure ||(a) shows the ground state condensation en- 
ergies for both even and odd grains calculated with 
g.c. BCS method, the PBCS approach (for N < 600) 
and Richardson's exact solution. The result E^ c M is 
also shown for comparison. The g.c. curves suggest the 
aforementioned "breakdown of superconductivity" ||||] 
at some critical p-dependent level spacing dp c above 
which E^ c = 0. In contrast, the £? BCS 's (i) are signif- 
icantly lower than the E^ s, thus the projection much 
improves the variational ansatz; and (ii) are negative for 
all d, which shows that the system can always gain en- 
ergy by allowing pairing correlations, even for arbitrarily 
large d. The exact solution £^ xact further improves the 
PBCS results, especially for intermediate level spacings. 
The PBCS results are evidently quite accurate for d A 
and, like the g.c. results, for d <C A. 

As anticipated in @ , the "breakdown of superconduc- 
tivity" is evidently not as complete in the canonical as in 
the g.c. case. Nevertheless, some remnant of it does sur- 
vive m 

£PBCS 

, since its behaviour, too, changes markedly 
at a p (and A) dependent characteristic level spacing dp 
(< dp* c ): it marks the end of bulk BCS-like behavior for 
d < dp , where Ep is extensive (~ and the start of 

a fluctuation-dominated plateau for d > dp , where Ep 
is practically intensive (almost d independent). 

The standard heuristic interpretation |14|] of the bulk 
BCS limit -A 2 /(2d) (which is indeed reached by E% 
for d — > 0) hinges on the scale A: the number of levels 
strongly affected by pairing is roughly A j d (those within 
A of sf), with an average energy gain per level of — A/2. 
To analogously interpret the d independence of Ep in the 
fluctuation-dominated regime, we argue that the scale A 
loses its significance - fluctuations affect all no = u>o/d 
unblocked levels within lod of Ef (this is made more pre- 
cise below) , and the energy gain per level is proportional 
to a renormalized coupling —Xd. 

The exact results smear out the crossover even more 
than the PBCS results (E™ act lacks the kinks of £^ BCS ), 



so much so that no sharply-defined crossover level spac- 
ing can be associated with £p Xact . However, the crossover 
scale evidently still is d ~ A. This can be confirmed by 
analyzing the functional dependence of the ground state 
energy on the coupling strength A: In the BCS limit, 
Ep — A 2 /(2g?), where A depends exponentially on A 
[since A = wr> sinh(l/A)]. In the fluctuation-dominated 
regime, however, perturbation theory in A suffices and 
the correlation energy is roughly linear in A. For each 
d, we thus fitted the numerical results for i?p Xact (d, A), 
calculated for various A, to 



E (X,d) 
E Q {X ,d) 



a(d) 



sinh(l/A ) 
sinh(l/A) 2 



13(d) 



(9) 



a "phenomenological ansatz" which intends to capture 
the relative importance of the exponential or linear A 
dependence in the coefficients a(d) and (3(d). The re- 
sults, shown in Fig. ||(b), clearly show the crossover from 
the BCS-dominated regime (a > j3) to the fluctuation- 
dominated regime (a < (3). 



D. Wave Functions 

Next we analyze the ground state wave- function, which 
can be characterized by 

C 2 Ad) = (cj. c j+ ct Cj -_) - (4 + c j+ )(ct Ci _), (10) 



a set of correlators that measure the amplitude enhance- 
ment for finding a pair instead of two uncorrelated elec- 
trons in a single-particle niveau \j, ±). For all j of an 
uncorrelated state one has Cj = 0. For the g.c. BCS case 
Cj = UjVj and the Cj's have a characteristic peak of 
width ~ A around ep, see Fig. |^(c), implying that pair- 
ing correlations are "localized in energy space" . For the 
BCS regime, both canonical methods produce Cj's virtu- 
ally identical to the g.c. case, vividly illustrating why the 
g.c. BCS approximation is so successful: not performing 
the canonical projection hardly affects the parameters Vj 
if d <C A, but tremendously simplifies their calculation. 
However, in the fluctuation-dominated regime d~ A, the 
character of the wave-function changes: weight is shifted 
into the tails far from sf at the expense of the vicin- 
ity of the Fermi energy. Thus pairing correlations be- 
come delocalized in energy space (as also found in [flBI ), 
so that referring to them as mere "fluctuations" is quite 
appropriate. In the extreme case d 3> A superconducting 
fluctuations are roughly equally strong for all interacting 
levels. 



E. Parity Effect 

The parity effect predicted in the g.c. ensemble 
can be studied with a simple generalization H of the 
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above methods. Specifically, we shall study the parity- 
dependent "pair-breaking energy" O p , i.e. the minimum 
energy required to break a pair by flipping a spin in an 
even or odd grain, defined as Qq = \{Ei — E ) an d 
fli = ^(E 3 — Ei), where E u denotes the energy of the 
lowest-lying state with u unpaired electrons with the 
same spin. (The pair-breaking energies can readily be 
measured in BRT's experiments by applying a magnetic 
field, whose Zeeman energy favors the breaking of pairs.) 
For a non-interacting system and approximatly also for 
large d ^> A, Hq ~ d/2 while f^i ~ d. On the other 
hand, in the bulk limit flo = fii = A, since in the bulk 
E p+2 = E p + 2 A. The parity effect now states that as d 
increases from d ~ 0, the pairing correlations die faster 
for an odd than an even grain, causing fii to initially de- 
crease faster than 51o. Since in the large d-limit f2i > flo, 
the two energies must cross somewhere at d ~ A, as 
shown in Fig. ||(d). This crossing is a manifestation of 
the parity effect already predicted in the g.c. framework. 

Remarkably, despite the crudeness and incorrect treat- 
ment of fluctuations of the g.c. method, it gives surpris- 
ingly good results for energy differences like the £l p . The 
fluctuations which it neglects seem to cancel substan- 
tially in energy differences such as fl p . This observation 
a posteriori justifies the use of the g.c. method even 
for small grains, at least for rough calculations of energy 
differences. 



IV. CONCLUSION 

In summary, the crossover from the bulk to the 
fluctuation-dominated regime can be captured in full us- 
ing a fixed- N projected BCS ansatz, or even exactly us- 
ing Richardson's method. With increasing d, the pair- 
ing correlations change from being strong and localized 
within A of £p, to being mere weak, energetically delo- 
calized "fluctuations" ; this causes the condensation en- 
ergy to change from being extensive to intensive (modulo 
small corrections). Thus, the qualitative difference be- 
tween "superconductivity" for d~A, and "fluctuations" 



for d~A, is that for the former but not the latter, adding 
more particles gives a different condensation energy; for 
superconductivity, as Anderson put it, "more is differ- 
ent" . 
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Note added: 

After submission of this manuscript works by Dukelsky and Sierra Jl(|[l7j used DMRG methods to numerically 
calculate the ground states of the presented Hamiltonian ([!]). Their results agree within a relative error of 1CP 4 with 
those obtained by Richardson's exact method. 
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FIG. 1. RBT's experimental tunneling spectra The distances between lines give the fixed- iV excitation spectra of the 
same grain containing (a) an even and (b) an odd number of electrons, as function of magnetic field. 
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FIG. 2. (a) The ground state correlation energies E p for even and odd systems (even: lower curves, odd: upper curves), 
calculated grand-canonically (GC), with PBCS and Richardson's method (exact) as functions of d/A — 2sinh(l/A)/iV. (b) 
shows the coefficients a and fi of Eq. (|!^) as discussed in the text. The pairing amplitudes Cj of Eq. ([To]) are shown in (c) 
for bulk (no symbols, dot-dashed), PBCS ("+", dashed) and exact results ("x", solid) for d = 1.09A and d = 4.34A. For 
d — 0.27A, the PBCS and exact curves (not shown) are indistinguishable from the bulk curve, (d), finally, shows the results 
the pair-breaking energies flo/i as calculated with g.c. BCS (dot-dashed) PBCS (dashed) and exactly (solid). (For clarity the 
graphs are offset by successively one unit each.) 
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